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Abstract 

The 1-Steiner tree problem, the problem of constructing a Steiner minimum tree 
containing at most one Steiner point, has been solved in the Euclidean plane by Geor- 
gakopoulos and Papadimitriou using plane subdivisions called oriented Dirichlet cell 
partitions. Their algorithm produces an optimal solution within Oiri^) time. In this 
paper we generalise their approach in order to solve the fc-Stciner tree problem, in which 
the Steiner minimum tree may contain up to k Steiner points for a given constant k. 
We also extend their approach further to encompass arbitrary normed planes, and to 
solve a much wider class of problems, including the fc-bottleneck Steiner tree problem 
and other generalised fc- Steiner tree problems. We show that, for any fixed fc, such 
problems can be solved in time. 

Keywords: fc-Steiner tree; bottleneck Steiner problem; network optimisation; polynomial- 
time algorithm 



1 Introduction 

The Steiner tree problem, which asks for a network with minimum total edge length inter- 
connecting a given set of points (called terminals), is a well known geometric variant of the 
spanning tree problem. The Steiner tree problem can be viewed as belonging to a family 
of problems where the aim is to construct a network T interconnecting a given set of n 
terminals, with the following properties: 

1. T may contain nodes other than the given terminals; 

2. T minimises a given cost function; 

3. the given cost function guarantees that T can be assumed to be a minimum spanning 
tree on its nodes (for a given metric). 



In addition to the Steiner tree problem under various metrics, this family also includes the 
power-p Steiner tree problem, where the cost of each edge of the network is the pth power 
of its length, and the bottleneck Steiner tree problem, where the cost of the network is the 
length of its longest edge and there is a bound on the number of extra nodes in the network. 
Both of these variants on the Steiner tree problem have numerous applications, particularly 
in the design of wireless communications networks. 
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Although the spanning tree problem is polynomially solvable, being solvable in O(n^) time 
in a general metric and in O(nlogn) time in the Euclidean plane, see [20], the problems 
in this more general family are mostly A'P-complete, even in the plane. This essentially 
stems from the fact that as the number of extra nodes that can be added to the network 
increases there is an exponential explosion in the number of different topologies that need 
to be considered. A natural way of controlling this increase in complexity is to bound the 
number of extra nodes, in other words, replace Property 1 above by the following: 

la. T may contain up to k nodes other than the given terminals, where A; is a given 
positive integer. 

We refer to problems in this modified class as generalised k-Steiner tree problems. Under 
suitable conditions, such problems are polynomially solvable. Understanding the efficiency 
with which such problems can be solved, when terminals are embedded in a normed plane, 
is the main topic of this paper. 

One of the seminal papers on this topic, for the 1-Steiner tree problem in the Euclidean 
plane, is a paper by Georgakopoulos and Papadimitriou [6], where an C'(n^)-time solution is 
given. They conclude their paper with a tantalising comment relating to their unsuccessful 
attempts at generalising their methodology to the A;-Steiner tree problem, even for k = 2. 
This comment has been a motivating factor for the current paper. Also, in their paper the 
proofs of some of the primary results are omitted, and some lack sufficient detail; another 
aim of this paper is to add some rigour to the more fundamental of these results. 

We may restate the goal of the Georgakopoulos and Papadimitriou paper as follows: find 
the point in the Euclidean plane which, if added to a given set of points, will result in the 
shortest possible spanning tree. The authors observed that one can significantly reduce the 
time-complexity of an algorithmic solution to the problem by first constructing a special 
partition. Given a set X oin terminals in the Euclidean plane, it is possible to partition the 
plane into 0{v?) regions such that if any new point s is embedded in the plane within one of 
these regions, say i?, then any minimum spanning tree T on XU {s} will have the following 
property: the neighbours of s in T will belong to some subset of a set Cx{R) containing at 
most six points from X, where Cx{R) is fixed for the given region R. This useful partition 
is referred to as the overlayed oriented Dirichlet cell partition, or OODC partition. Their 
algorithm takes as input the set X of terminals and then starts by calculating the OODC 
partition, the set Cx{R) for every R, and a minimum spanning tree T' on X. All this, as 
well as a preprocessing step on T', is done within a time of 0{n'^). For each region i?, the 
algorithm then iterates through all subsets S of Cx{R) and calculates s, the Steiner point 
of the nodes in S (this takes constant time for each S). The algorithm then updates T' 
(which can also be done in constant time because of the preprocessing step) to include s. 
The cheapest tree is selected at the end as an optimal solution. A naive algorithm for the 
1-Steiner tree problem would attain a complexity of 0(n®.n log n) since it would have to 
iterate through all subsets of up to six terminals and then calculate the optimal position of 
the Steiner point and the corresponding minimum spanning tree for each of these subsets. 
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The powerful simplifying properties of the OODC would clearly be advantageous in a 
generalisation of the algorithm to arbitrary normed planes. However, the construction of 
the OODC partition given in [BJ is valid for the Euclidean plane only. We will provide 
a new method, based on abstract Voronoi diagrams, for constructing the partition for 
terminals embedded in an arbitrary normed plane. The construction is based on theoretical 
results, but we also define a class of norms for which the algorithm would be practically 
implement able. Once the partition has been found, our algorithm calculates the optimal 
positions of all k Steiner points simultaneously. Of course, these positions will depend not 
only on the neighbours of the Steiner points, but also on the cost function of the given 
generalised Steiner tree problem. Since, at the start of this step, the neighbours of the 
Steiner points are fixed but the Steiner points are free, this subproblem is a generalised 
version of the well-known fixed topology Steiner tree problem (discussed in Section H]). Even 
though we consider k to be part of the input of our algorithm, k will generally be much 
smaller than n, and therefore we assume that this step can be done in constant time. A 
novel method for updating a minimum spanning tree is then utilised to calculate a potential 
solution for every choice of coordinates of the Steiner points. Once again, a cheapest tree 
is selected as the optimal solution. The total time-complexity turns out to be 0{n'^^) when 
factors that are not functions of n are excluded. 

There are a number of authors who have looked at adapting the solution to the 1-Steiner 
tree problem in [6] to other ip norms. Kahng and Robins in [12] do this for the rectilinear 
plane, however, their paper only uses the solution as a step in a heuristic algorithm for the 
rectilinear Steiner tree problem, and not much attention is devoted to the solution of the 
1-Steiner tree problem itself. Griffith et al. expand on this heuristic idea in the rectilinear 
plane, see 0, and they provide a simple procedure (though without proof) for updating 
a minimum spanning tree when a new node is introduced. Lin et al. in turn adapt the 
approach presented by Kahn and Robins to the Aa-plane in |15j . 

Section [2] provides some preliminary definitions. We define the OODC partition in Section 
[3] and provide details, including the time-complexity, of its construction. We also present 
three restrictions on the given normed plane that would allow the construction to be im- 
plemented in practice. The problem of locating Steiner points under a fixed topology is 
outlined in Section HI Many results on this topic exist in the literature, so we summarise 
them very briefly for a few of the more common instances of the generalised Steiner tree 
problem. Section [5] describes a new method, and furnishes a correctness proof, of updating 
a minimum spanning tree to contain a fixed subtree. In Section [6] we present our algorithm 
for calculating an optimal solution to any given instance of a generalised A;-Steiner tree 
problem in a normed plane, and then prove its correctness and verify its time-complexity. 

2 Preliminaries 

We begin by formalising the definition of a generalised /c-Steiner tree problem, sketched in 
the introduction. 

Throughout this paper we use the symbols E{G) and V{G) for the edge-set and node-set 
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respectively of a graph G. We also use the notation G = {V{G),E{G)). Let A;' > be 
given; let || • || be a given norm on M^. 

Given a set P = {pi, . . . ,Pn+k'}i let {T} represent the set of all spanning trees for the 
elements of P. For each T there is a corresponding set of edges E{T) = {ei, . . . , Cn+k'-i} 
(with each Ci G P x P). Let X = {xi, . . . , with each Xi eM?, represent an embedding 
of the set {pi, ■ ■ ■ ,Pn} in (where each Xi is an embedding of the corresponding pi). 
We can think of each T as representing the topology of a tree network interconnecting 
X and using k' extra nodes, and we can equate the edges E{T) with the arcs of such a 
network. For a fixed embedding of this network we let S = {xn+i, ■ ■ ■ ,Xn+k'}-i with each 
G M^, be the locations of the extra nodes corresponding to {pn+i, ■ ■ ■ ,Pn+k'}- We refer 
to X as the set of terminals and 5 as the set of Steiner points of the network. Now let 
^T,x,s = (Ikill, • • • , ||e„_|_/;/_i||); i.e., the components of gt,x,s are the edge lengths of such 
a network, for a given tree topology and a given set of embedded nodes. Such a vector is 
well defined up to the order of its components. 

Let a : W]^^ "''^ — >■ R be a symmetric function (ie, independent of the order of the compo- 
nents of the vector on which it acts). We think of ck as a cost function on a tree network. 

In other words, a{G'Y^x,s) is the cost of the network with topology T and nodes X and S*, 
and vam.a{e-Y^x,s) is the minimum cost (with respect to a) of any tree interconnecting the 

nodes X and k' other points. Hence, for the power-p Steiner tree problem we define 

n+k'-l 

a{er,x,s) = ap{^T,x,s) ■■= XI H^^H^' 

1 

whereas, for the bottleneck problem (where the cost of the network is the cost of the longest 
edge) we have 

a{er,x,s) = aoo(er,x,s) = max [|ei||. 

E{Tj 

We define such a symmetric function a to be ii-optimisable if and only if there exist T* 
and S* such that a{e-j-* ^x,s*) = ^^''^Oi{^T,x,s) and ceilej-* ^x,s*) = niinai(e7-^x,s*)- In 

other words, a is £i-optimisable if for any given X there exists a tree T interconnecting X, 
with minimum cost with respect to a, that is a minimum spanning tree on its complete set 
of nodes. We denote the cost of such a tree by HTHq,. It is easy to show that ccp, for p > 0, 
and a^o are ^i-optimisable. 

Definition. For any given positive integer k, a generalised k-Steiner tree problem is defined 
to be any problem of the following form: 

Given A set X o{ n points in M^, a norm || • ||, and a symmetric £i-optimisable 
function a. 

Find A set S of k' < k points in R^, and a spanning tree T on X Li S with 
topology T such that ||r||a = a{er,x,s) = min a(er',x,5')- 
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We refer to T as a generalised k-Steiner minimum tree. The next lemma is an extension of 
the Swapping Algorithm, found in [14]. It follows from the matroid properties of minimum 
spanning trees. 

Lemma 1 Let T be a minimum spanning tree on the terminal set X, and let T' be a 
spanning tree for X. We can transform T' to T by a series of edge swaps, where each swap 
involves replacing an edge Cj G E[T') by Cj € E[T) such that ||ej|| > ||ej||. 

The corollary below shows that T is equivalent to any minimum spanning tree on X U S". 

Corollary 2 If T' is a generalised k-Steiner minimum tree on X with Steiner points S 
then every minimum spanning tree on X U S is a generalised k-Steiner minimum tree on 
X. 

Proof. Let T be a minimum spanning tree on X U 5*. By Lemma [H we can transform T' 
to T by a series of edge swaps, each of which replaces an edge with another of the same 
length. By the symmetry of a each such edge swap does not increase ||T'||q,. ■ 

Throughout this paper we perform various constructions involving the unit ball B for the 
given norm || • ||, for instance calculating the intersections of two unit balls. Our main 
interest in this paper is in the computational nature, specifically the time-complexity, of a 
solution to any instance of the generalised /c-Steiner tree problem. Li order to find efficient 
algorithms, we need to compute these unit ball operations to within any fixed precision in 
constant time. We therefore restrict the norm || • || so that its unit ball is always simple 
enough to perform these operations. We will provide more detail regarding these restrictions 
in the next section. For similar computational reasons we will also be placing a restriction 
on a, and this is discussed in Section [H 

3 The Overlayed Oriented Dirichlet Cell Partition 

Let a norm 1| • || on be given with corresponding unit ball B. Our aim in this section 
is to describe the construction of the oriented Dirichlet cell (ODC) partition for any set 
X oi n terminals embedded in this normed plane, and to show that it can be constructed 
within a time of O(nlogn). We also show that, with a time complexity of ©(n^), multiple 
ODC partitions can be overlayed. This final partition is the afore-mentioned overlayed ODC 
partition (OODC partition), and is a core component of our algorithm. 

Georgakopoulos and Papadimitriou [6] allude to a simple method of constructing an ODC 
partition for terminals embedded in the Euclidean plane. Unfortunately this method does 
not work for arbitrary normed planes. We circumvent this problem by defining a type of 
abstract Voronoi diagram that is equivalent to the ODC partition, and then showing that 
this Voronoi diagram can be calculated in the required time. 
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We now state the first of three restrictions on B. We defer a discussion of these restrictions 
(including the description of a class of norms that satisfy all of them) to the end of the 
section. Let the boundary of B be denoted by hd{B). 

Restriction 1 The intersection points of any two translated copies of bd(i?), and the 
intersection points of any straight line and hd{B), can be calculated to within any fixed 
precision in constant time. 

Lemma 3 There exist six points {yi : i = 0, ...,5} on hd{B) such that for any pair of 
rotationally consecutive ones, say yi,yj, we have \ \yi — yj|| = 1- Moreover, these six points 
are constructive . 

Proof. The standard ruler and compass construction of the hexagon will produce these 
six points, where B plays the role of the circle in the construction. Given any point ys on 
bd(-B) construct a translation of bd(i?) centered around ys. Let yo be the first intersection 
point of the two boundaries as we traverse the boundary of the original ball anticlockwise 
from ys. Let yi = yo — Ys- Note that this point also lies on bd(i3) and that yoyioys is a 
parallelogram. The remaining three points are constructed using the central symmetry of 
B. See Figure [T] for an example where i? is a tilted ellipse. The distance properties of the 
lemma follow easily by construction. ■ 




Figure 1: The standard hexagon construction 

For any two directions c^j and (t)j in the plane, we denote by K{y,(j)i,(j)j) the cone with 
vertex y from the ray with vertex y in direction (pi to the ray with vertex y in direction <j)j. 
For each y^ from Lemma [3] let 9i be the direction of the ray oy^, where o is the center of 
B. We assume that the {6i} are ordered in an anti-clockwise manner, and two consecutive 
directions will be denoted by 9i and 6i+i (i.e. the mod 6 notation will be omitted). As 
another example we show, in Figure [21 the six directions produced when B is the unit ball 
of the rectilinear plane. 

Lemma 4 Let x ^ B, x ^ o and a and b points on the boundary of B such that the 
segments ao and bx intersect in P. Then ||a — x|| < ||b — x||. 
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Figure 2: The unit ball of the rectilinear plane and corresponding {9i} 

Proof. Applying the Triangle Inequality to Apax and Apbo, we obtain ||a — x| | + ||b|| < 
||b — x|| + ||a|| from which the lemma follows; see Figure [3l ■ 



b 




Figure 3: Illustration of the proof of Lemma H] 

Lemma 5 Let y be any point in the plane. Then there exists a minimum spanning tree T 
on X [J {y} with the following property: for each i = 0, 5 there is at most one point of 
X adjacent toy in T and lying within K{y,6i,6i^i), and this point is a closest terminal to 
y in the cone. 

Proof. Let T' be a minimum spanning tree on X U {y}. Let xq G X be a terminal in 
K{y,6i,6i+i) that is closest to y, and suppose that (y,xi) G E{T') where xi G X is any 
other terminal in K(y,6i,6i^i). We show that we can replace the edge (y,xi) in T' by 
either (y,xo) or (xi,xo) so that the resulting tree is still a minimum spanning tree on 
X U {y}. From this, the statement of the lemma follows. 

Assume first that the path in T' connecting y and xq passes through xi. In this case we 
can replace (y,xi) by (y,xo) without losing connectivity or increasing the length of T'. 
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Assume, on the other hand, that the path in T' connecting y and xq does not pass through 

Xi. 

Claim: ||xo — xi|| < ||y — xi||. 
Without loss of generahty, we assume that y = o, l|xi|l = 1, and K{o,9i,9i^i) intersects 
the boundary of the unit ball B in an arc from a to b (with ||a — b|| = 1). We can also 
assume, without loss of generality, that xi lies on the same side of the line through oxq 
as b. The convexity of B implies that the line segments oxi and ab intersect, hence by 
Lemma |4] we have 

||a - xill < l|a - bll = 1. (1) 
We now prove the claim via two cases, illustrated in Figure HI Firstly, suppose xq and o 




Figure 4: The two cases of the Claim in the proof of Lemma [5] 

are on the same side of axi (including the case where xq lies on axi). By Inequality ([1]) a 
lies in the unit ball centered at xi, so, by convexity, xq also lies in this unit ball. Hence, 
||xo — xi[| < 1 as required. For the second case, suppose that xq and o are on opposite 
sides of axi. Let the ray from xi passing through xq intersect B at x.'q. Then Xq and o are 
also on opposite sides of axi, and hence, by Lemma [H ||xo — xi|| < ||xq — xi|| < ||a — xi||. 
Therefore, ||xo — xi|| < 1 by Inequality ([T|), and the claim is proven. 

By the claim we can now replace the edge (y,xi) by (xi,xo) without losing connectivity 
or increasing the length of T'. ■ 

For each i = 0, 5, the ith oriented Dirichlet cell (ODC) of w G X is the set: 

{y e M2 : ||w - y[| = min{||x - y||,x G X n E:(y, 0^, 

In other words, this is the set of all points {y} whose closest terminal in the cone K{y, 6i,6i^i) 
is w. We will show that the set of ith ODCs, called the ith ODC partition of X is a type 
of Voronoi diagram. 

In [2] Chew and Drysdale present an "expanding waves" view of Voronoi diagrams. If n 
pebbles are dropped simultaneously into a pond, the places where wave fronts meet define 
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the Voronoi diagram on the n points of impact. In the EucUdean case the wavefronts are 
circular, but in theory any closed convex oriented curve can qualify as a wavefront and 
thereby define an abstract Voronoi diagram. For any such shape C and set of terminals X 
we say that the resulting diagram is the Voronoi diagram of X based on C. 

We may define this Voronoi diagram more formally as follows. Let (5c : ^ be the 
distance function based on C; in other words, for any points xq, xi € we let ^cC^O) ^i) = 
inflA""*^ : A(xi — xq) € C}. We then define a region 14; = {y : (5c(x, y) G M and 5c(x, y) = 
min{5c(x',y) : x' G X}^ for each x G X. The set {Vy^} is the required Voronoi diagram 
based on C. 

Proposition 6 For any i = 0, ..,5 the ith ODC partition of X is equal to the Voronoi 
diagram of X based on the sector B n ^(o, 180° + 9i, 180° + 6*4+1). 

Proof. This follows immediately from the central symmetry of i3. ■ 

In Figure [5] we give an example of an ODC partition when the original unit ball is a circle 
(i.e. the Euclidean case) and therefore C is a circular sector, and in Figure [6] we give an 
example of what the boundary of wavefronts look like when C is a regular hexagon. 
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Figure 5: An ODC partition of a three-terminal example 
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Figure 6: Regular hexagon based Voronoi diagram for two points 
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The next theorem now gives us the required time-complexity for calculating the ith ODC 
partition under certain conditions. 

Theorem 7 The Voronoi diagram of n points based on a closed convex shape C can he 
constructed in 0{n\ogn) time and 0{n) space as long as the following operations can he 
performed in constant time: 

1. Given two points, find the boundary where the two wavefronts meet. 

2. Given two such boundaries, compute their intersection (s). 

We therefore also impose the following restriction on B. 

Restriction 2 Let C be any sector of B. Then, given any two points, we can find the 
boundary where the two wavefronts based on C meet, and, given two such boundaries, we 
can compute their intersection. Moreover, these operations can he performed to within any 
fixed precision in constant time. 

The next step is to overlay the six ith ODC partitions. The theorem we use, which is 
a result from [6j, assumes that the regions in each partition have boundaries consisting of 
straight line segments. We therefore state our third and final restriction on B. 

Restriction 3 The shape of B implies that the ith ODC partition of any set of points is 
linear (i.e. the boundary of any ODC consists of straight line segments). 

Theorem 8 Let q be any positive integer. Then q linear plane partitions can be over- 
layed in 0{q^n'^) time, where n is the total number of regions in each partition. 

As a consequence of the previous results, within O(nlogn) total time we can calculate the 
ith ODC partition of the plane for each i, and then, in a time of O(n^), overlay these six 
partitions resulting in the OODC partition. It is easily observed that the OODC partition 
has O(n^) regions. 

Let i? be a region of the OODC partition. Note that R is the intersection of at most 
six ODCs. Let {Dj} be the set of these ODCs, where Dj is an ODC of pj € X. Let 
Cx{Fi) = {pj}. The power of the OODC partition lies in the next theorem, which now 
follows from Lemma [5l 

Theorem 9 Let s be any point in R. Then there exists a minimum spanning tree T on 
X U {s} such that the set of neighbours of s in T is a subset of Cx{R)- 

The question arises as to whether norms exist with unit balls satisfying all three restric- 
tions. Let V be the class of norms defined by the condition that each norm's unit ball is 
either a polygon or an ellipse. Suppose that || - || G V and that the corresponding unit 
ball is B. Clearly Restriction 1 is true for B. Given any two points, their bisector (based 
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on a sector of B) will be a polygonal line that can be computed with some simple vector 
operations (see [2]). The same holds true for intersecting two boundaries, and therefore 
Restriction 2 is satisfied. Since any ith ODC partition consists of segments of bisectors and 
segments of the limiting rays of i^(y, Oi, ^j+i) for some y, Restriction 3 follows immediately. 

4 Generalised Steiner Tree Problems for a Fixed Topology 

Each main iteration of our algorithm produces a set of neighbours (from X Li S) for each 
of the k' < k Steiner points. Finding the optimal coordinates of the Steiner points for 
the subtree induced by the Steiner points and their neighbours is a problem known in the 
literature as the fixed topology Steiner tree problem. We state the problem more formally 
as follows: given a set X' of c < Qk' embedded terminals, a set S of k' free (i.e. non- 
embedded) Steiner points, and a tree topology T spanning all these nodes, we wish to find 
the coordinates of the Steiner points (i.e. find the set S) such that a{e'T-^x' ,s) is minimised, 
where Gf x'.s = (l|ei||, ||ec+fc'-i||)- For our purposes a node x oiT will be of degree one 
if and only \i x ^ X' . This problem is interesting in its own right, but is also a key step of 
our main algorithm. We therefore restrict a so that this problem is always computable to 
within any fixed precision in constant time, for a fixed value of k. In fact, k (and therefore 
c) will generally be much smaller than n, and therefore we are not particularly interested 
in the time-complexity of this step. 

We discuss briefly a few functions that satisfy this restriction. 

(1) a{e'j-^x',s) = X] I Nil I- t'^is case we are dealing with the well-known Steiner tree 
problem for a fixed topology. In the Euclidean plane the problem has an C'(c^)-time 
solution provided that no point has degree larger than 3, see [10]. Unfortunately, for 
the /c-Steiner tree problem degree 4 points do exist (but degree 5 do not; see [18]). 
and this more general problem may be NP-complete. A similar result holds for the 
rectilinear and other fixed orientation planes [l]. 

(2) a{eq-^x',s) = ll^ill'^i p > 0. This is referred to as the power-p Steiner tree problem 
for a fixed topology. In the Euclidean plane with p = 2, Ganley [1] shows that the 
problem can be solved within time 0{c). 

(3) a{eq- X' s) = Ihii fy^||ej||^^ , i.e. the bottleneck Steiner problem for a fixed 

topology. This problem has an 0{(?) solution in the rectilinear plane, see [5j. In the 
Euclidean and general ip planes there exists various numerical algorithms that can 
calculate a solution to any desired precision, see for instance [3], [16]. A fully poly- 
nomial time approximation scheme (FPTAS) exists for the problem in the Euclidean 
plane (see [4]). 
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5 Minimum F-fixed Spanning Trees 



Many papers exist in the literature that deal with the time-complexity of updating a mini- 
mum spanning tree when a new node is introduced; see for instance [llj where they show 
that a tree on n nodes can be updated with a new node in 0(logn) parallel time using 
n/logn exclusive read, exclusive write, parallel random access machines (EREW PRAMs). 
Georgakopoulos and Papadimitriou utilise a preprocessing step in ^ so that a minimum 
spanning tree can be updated in constant time with a new point. 

In order for our algorithm to find an optimal tree for all Steiner points simultaneously, it 
must update a minimum spanning tree (in constant time) to include a new subforest F. 
However, it is not necessary for the updated tree Tp to be a minimum spanning tree on 
its nodes, but only that it is a shortest total length tree spanning X U V{F) such that the 
neighbour-set of each Steiner point in F remains unchanged in Tp. 

Let N{T, s) denote the set of neighbours of a node s in a graph T. A forest F with 
node-set X' U S, where X' C X and S" C with 15"! < k, is called feasible if and only if 
{x G V{F) : X is a leaf of F} = X' and \N{F, s)| < 6 for every s G 5. If 5 is the set of at 
most k new embedded Steiner points then, in our main algorithm, the forest F induced by 
the set of new edges incident to elements of S will be feasible. A shortest total length tree 
Tp, such that V{Tf) = X \J S and N{Tp,s) = N{F,s) for every s G S, is referred to as a 
minimum F-fixed spanning tree. 

We use the symbol y) to represent a path through T with endpoints x and y, and 

we use y) to denote the longest edge on Pr(x, y). We will make use of the following 
theorem. 

Theorem 10 !ll3f A tree T is a minimum spanning tree on X if and only if for every pair 
of non- adjacent nodes x,y G X, \\e\\ < ||x — y|| for every e G E{PT{x,'y)). 

Let T now be a minimum spanning tree on X. The preprocessing stage PSl calculates 
^^(x, y) for every pair of nodes x, y G V{T). This requires O(n^) time and 0{n?) space. 
We incorporate a consistent tie-breaking component for choosing between edges of exactly 
the same length by placing an order on E(T) and choosing the earlier edge in this ordering 
whenever a tie occurs. 

Proposition 11 Let T be a minimum spanning tree on X, and assume that the above 
preprocessing step has been performed. If F is connected and feasible, then a minimum 
F-fixed spanning tree Tp can be constructed from T in 0{k'^) time. 

Proof. Let G = T U F, let A = V{F) n X, and note that \A\ < 6k. A number of cycles 
may occur in G, each one of them containing a path through F with endpoints from A. Let 
T' be the graph obtained by deleting the set of edges {£T(xi,Xj) : Xi,Xj £ A,i j} from 
G. We will show that T' = Tp, which suffices to prove the proposition since T' can clearly 
be constructed in 0{k'^) time. 
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To prove that T' = Tp we first show that T' is acychc and spans X U 5. Every cycle 
of G is of the form Pp{xi,Xj), PT{xj,Xi), and therefore deleting every £r(x,j,Xj) from G 
produces an acyclic graph. We use induction on |^| to prove that T' is connected. Let 
Ah = {xi, .., Xf,} C A for some b S {2, .., 6k}, let Fh be the subtree of F induced by 5 U A^, 
and let G^ = T U Ff,. Subtracting L5 = {£T(xi,Xj) : 1 < i < j < b} from E{Gh) produces 
the graph T^. For the base case we let 6 = 2. The only cycle of G2 is Pf(xi, X2), Pt(x2, xi), 
and ^t(x2,xi) is an edge of this cycle. Therefore deleting £t{^2,^i) does not destroy the 
connectivity of T2 on X U 5. 

Next assume that Tj, spans X U 5 for some 2 < b < Qk — 1 and suppose that Xf,+i G X\Ai,. 
Since Tb is connected and acyclic there is exactly one path connecting Xf,+i to a node of Ah 
not passing through any element of S, i.e. this path is of the form PT(xb+i,Xj.) for some 
unique x^ € A^. Let A^+i = U {x5_|_i} and let s = A^(F, x^+i). Then T{,+i is the graph 
with V{Th+i) =XUS and E{n+i) = {E{n) U {(s,X6+i)}) \{^T(xb+i, x^) : x, G A}- 

Claim: For every x, G Ai, either ^^-(xfc+i, Xj) G or £T(x6+i,Xj) = ^T(xfc+i,Xr). 
Let Xj G j4fe\{xr} and consider the following two cases. If x^+i lies on P7^(xj,Xr) then 
£T(xi,Xr) = max{£T(xi,Xb+i),^T(xb+i,Xr)} = £T(xi,Xfc+i) since Pr(xb+i,Xr) is a path in 
Tb and therefore does not contain ^^(xjjXr). Therefore £T(xf,+i,Xj) G Lb- For the second 
case, if x^+i does not lie on Pxi'^i^'i^r) then let y be the first common point of the paths 
PT(xi,Xr) and PT(x6+i,Xr); see Figure El 




Figure 7: The second case of the claim 

Note that y may be equal to x,.. Clearly 

^T(xi,y) >^T(y,x^) (2) 

since PT'(y,x,.) is also a path in Tb- There are now two possibilities to consider; either 
^T(x;,+i,Xi) = £T(xi,y) = £T(xi,Xr) G Lb, or iri^b+i^^i) = ■^T(x6+i,y) = £T(x6+i,Xr), 
where for each possibility the second equality follows from Inequality ([2]). The claim follows. 

By the above claim E{TbJ^i) = {E{Tb) U {(s, X5_|_i)}) \{^t(x;,_(_i, x^)} and we deduce that 
Tfc+i has been constructed from Tb by adding one edge from F and then deleting an edge of 
T on the resultant cycle. This completes the induction argument, and hence T' is connected 
and spans X U S. 
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Next we prove that T' is a minimum F- fixed spanning tree. Let K be the complete graph 
on X. Furthermore, suppose that the edges of K are weighted by the function w, where 



if X G A and y € ^, 

llx — vll otherwise. 



Let Ta be any spanning tree on A. Then clearly T' is a minimum fixed spanning tree 
if and only if the graph Tk, where V{Tk) = X and E{Tk) = {E{T') n (X x X)) U Ta, is 
a minimum spanning tree of K with the above weight function. But this follows from a 
simple application of Theorem I lOi Hence T' = Tp, as required. ■ 

An immediate consequence of the above proof is the following result. 
Corollary 12 {{iri^i^^j) ■ G ^}| = !^| - 1- 

Next we generalise the Proposition [11] to the case where F is not necessarily connected, 
which must be considered if k > 1. If A; > 1 we perform an additional preprocessing stage, 
PS2, to calculate a TRUE/FALSE table H such that i?e,y,z = TRUE if and only if edge 
e G E{PT{y,z)). This requires at most O(n^) time and 0{n'^) space. For each connected 
component of F let A" = V{F'-) n X. We claim that Algorithm F-MST, given in Table 
[U calculates a minimum F-fixed spanning tree for any feasible forest F. 

In essence. Algorithm F-MST deletes the longest edge not in F, of every cycle in T U F. 
The shortest path through J connecting some pair x, y € A^ is simply a representation 
of the sequence of paths through T that are contained in Prpi-i{x,y). The longest edge 
of all paths in this sequence is therefore the longest edge on the path Fpi-i(x, y) barring 
any edges of F. It should be clear that any inductive method that correctly deletes the 
longest edge of each newly created cycle, is equivalent to deleting all longest cycle-edges 
simultaneously. To formally prove correctness of Algorithm F-MST we need the following 
lemma. 

i-l 

Lemma 13 Let x, y G A^. Then ^*(x, y) is the longest edge in F(Fri-i (x, y))\ |^ E{F^). 

i=i 

Proof. Let J be the graph defined in Step la(i) of Algorithm F-MST, and let P' = 
{x}, A^^ , ...,A^'^, {y} be any path in J, where clearly 1 <bj < i — 1. Note first that o"i and 
(72 are well defined functions on F( J ) since, as we will show in the next proposition, r*~^ 
is a tree. If ei, Cd+i is the edge-sequence of P' then there is a corresponding x — y walk 
through r*~^ of the form 

W(P = (^71 (ei ) , ^2 (ei ) ) , , Pt (ai (62) , ^2 (62) ) , . . . , Ft (^Ti (e^+i ) , ^2 (cd+i ) ) 

where: <ti(6i) = x; 0-2(6^+1) = y; (7i{ej) G A^^-^ ioi 1 < j < d + 1; cr2{ej) G A''^ for 
1 < j < d; every F ^ is a path in F^J between two of its leaves, or is the empty set. 



14 



Algorithm F-MST 



Input: A set X of points in the plane, a minimum spanning tree 
T on X, and a feasible forest F with t connected 
components. 



Output: A minimum i^'-fixed spanning tree Tp on X U V{F). 



Step 



Description 



1 

la 

la(i) 



la(ii) 
la(iii) 



Let £)0 = and let = T. 

For i = 1 to t Do 
Begin 

For every distinct pair x, y G Do 
Begin 

Let J be the graph with 

{{x},{y}} ifz = l, 

{{x},{y},^i,...,^-i} ifz>l. 

and 

E{J) = {{U, U'):3w eU A 3w' G U' such that 



V{J) 



= FALSE Ve G D''^}. 



For every e = {U,U') ^ E{ J) let (Ti(e) = w and let 
C72(e) = w' (where w, w' are from the previous step). 
Perform a search through J to find the path 
P = P/dx}, {y}) such that |£^(P)| is a minimum. 
Let ^*(x,y) be the edge from {£T(cTi(e), a2(e)) : e G E{P)} 
of maximum length. 
End 

Let L^ = {r(x,y) :x,y gA^}, 
D« = D^-i u L\ 

= {V{G'),E{G')\D'). 

End 

Let Tf = TK 



Table 1: Algorithm F-MST 
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Let W' be any x — y walk through T*"^, and suppose that as we trace W' from x to 
y, we consecutively intersect Steiner points of the distinct trees F^^ , F""" , where each 
I < aj < i—l. Let W~^(W') = {x}, A'^^ , A'^'' , {y}. Note that it is not necessarily the case 
that W-^'iWiP') = P'. Now suppose that P = Ty-^^T-i (x, y)) = {x}, ^'^^ , ^-^^ , {y} 
(which, note, is a path through J). Since W{P') contains the unique path Pj'i-i(x, y), the 
node-sequence of P must be a subsequence of the node-sequence of P'. Since this holds 
for any P', this implies that P is the path through J of minimum length (i.e. minimising 
|i?(P)|) connecting {x} and {y}. The lemma follows. ■ 

Proposition 14 Let T he a minimum spanning tree on X, and assume that preprocessing 
steps PSl and PS2 have been performed. If F is a feasible forest then Algorithm F-MST 
correctly produces a minimum F-fixed spanning tree Tp in 0{k'^^^^k\) time. 

Proof. The proposition is easily verified by using induction on t (the number of connected 
components of F). Proposition [11] proves the base case t = 1, and similar reasoning is 
used to prove that each subsequent T* is connected. At each inductive step, Lemma \T3\ 
assures us that Algorithm F-MST correctly deletes the longest edge (excluding edges of F) 
of any new cycle formed. To prove minimality of Tp we once again construct (analogously 
to Proposition [TT]) a weighted complete graph K on X and a tree Tk , such that Tp is a 
minimum F-fixed spanning tree on X U V{F) if and only if Tk is a minimum spanning tree 
of K. Theorem 1101 then completes the minimality proof. 

To verify the time-complexity first note that t < k and \ A^\ < 6k. Step la of the algorithm 
requires 0{k^) time. Step la(i) requires 0{k'^^k\) time and Step la(iii) requires 0{k) time. 
■ 

Recall that when /c = 1 we do not need to calculate the table H during preprocessing 
because Proposition [11] will apply. Therefore the space complexity reduces to 0{n'^) in this 
case. 

6 The Main Algorithm 

We present Algorithm k-GSMT, in Table [2] 

By using Cayley's formula and the observation that each spanning forest is a subgraph 
of a spanning tree which has k — 1 edges, we get an upper bound for fi{k) of 126^. /c'^"^ 
in Step 4b of the algorithm. The function /2 in Step 4b(i) will depend on the relevant 
generalised fc-Steiner tree problem. We mention once again that we consider k to be much 
smaller than n, and therefore we conclude that the overall time-complexity of Algorithm 
fc-GSMT is 0(n2*=). 

We now demonstrate the correctness of the algorithm. Let Topt be a generalised /c-Steiner 
minimum tree on X. Let S be the set of Steiner points in Topt and let Fopt be the sub-forest 
of Topt induced by the edges of Topt incident with elements of S. Let Top^ be a connected 
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Algorithm A;-GSMT 


Input: A set X of n points in the plane, a unit ball B, a positive integer k, 
and a symmetric ^i-optimisable function a. 


Output: A set S of at most k Steiner points, and a tree T* interconnecting 

X U S, such that \\T*\\a = mma{er,x,S')- 

T,S' 


Step 


Description 


Time 


1 
2 
3 


Construct the OODC partition of X. 
Construct a minimum spanning tree T on X. 
Perform preprocessing steps PSl and PS2 on T. 


0{n logn) 
O(nlogn) 
0{r?) 


4 

4a 


For every k' < k and each choice (with repetition) 
of k' regions, Rk', of the OODC partition Do 
Begin 

Associate the free Steiner point Sj with region Ri. 
Let Q be the graph consisting of the vertices 
[jCxiRi) U {si,...,Sfc/}, aU edges {si,Sj),i / j, and 
all edges (si,x) for every x G Cx{Ri)- 
Let Q* be the set of all feasible subforests of Q. 




4b 
4b(i) 
4b (ii) 

5 


For each J" G ^* Do 
Begin 

Solve the fixed topology generalised Steiner 
tree problem for T to get the tree F. 
Run Algorithm F-MST with input T and F, 
and let Tp be its output. 
End 

End 

Select a smallest total cost Tp produced and let T* = Tp. 
Let S be the set of Steiner points of T*. 


oum) 

OU2{k)) 

Oik'^'^+^kl) 



Table 2: Algorithm A;-GSMT 
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component of -Fopt with ki Steiner points and terminal set C X. Note that, hke Topt, 
-Fopt may not be unique for a given set X. We have the following two easy lemmas. 

Lemma 15 F^^^ is a generalised ki-Steiner minimum tree on A^. 

Lemma 16 Let be any generalised ki-Steiner minimum tree on Al . Suppose we trans- 
form Topt to T' by replacing the subtree F^^^ on T by Y^. Then T' is also a generalised 
k-Steiner minimum tree on X. 

Proposition 17 Algorithm k-GSMT constructs a tree T* that is a generalised k-Steiner 
minimum tree on the terminal set X . 

Proof. By the properties of the OODC partition, during the course of the algorithm a forest 
F induced by edges incident with Steiner points is constructed with connected components 
F^ such that each F^ has the same terminal set (i.e. A^) and the same topology as i^opt, 
and therefore 

11-^ \\a = ll-P'optlU- (3) 

We now consider two cases. 

Suppose, for the first case, that Fopt = F- Step 4b(ii) of the algorithm constructs a 
minimum F-fixed spanning tree Tp on X L) S. Since Topt is a minimum spanning tree on 
XUS and contains T as a sub-forest, it follows that Tp is also a minimum spanning tree on 
X L) S. By Lemma El T* = Tp is a generalised fc-Steiner minimum tree on X, as required. 

If, on the other hand, Topt 7^ T, then there is a tree Tp constructed in Step 4b(ii) of the 
algorithm that is the same as in the previous paragraph, except each T* is replaced by T^pt. 
By Equation ([3]) and Lemma [16] it again follows that T* = Tp is a generalised fc-Steiner 
minimum tree on X. ■ 

7 Conclusion 

The outcome of this paper is a generalisation, on multiple fronts, of Georgakopoulos and 
Papadimitriou's 0(n^) solution to the 1-Steiner tree problem. By utilising abstract Voronoi 
diagrams, we build on their complexity-reducing concept of oriented Dirichlet cell partitions. 
The result is a broadening of the scope of these partitions to include terminal sets in 
arbitrary normed planes. A bigger challenge in our research was to construct a generalisation 
to k Steiner points. We achieve this by producing a novel method of updating a minimum 
spanning tree to include a fixed subtree. A two-part preprocessing stage allows this to 
be done in constant time with respect to the total number of terminals. One of the key 
observations of our research was that the main algorithm from [6] basically pertains to any 
"Steiner- like" problem with cost function a, as long as it is guaranteed that a solution exists 
which is optimal with respect to a and is also a minimum spanning tree on its complete set of 
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nodes. This fact allows us to accommodate the class of generalised fc-Steiner tree problems 
with symmetric ^i-optimisable cost functions. The result is an C'(n^'^)-time solution to this 
class of problems. 

It may be possible to generalise our algorithm to higher dimensional spaces, at least in the 
Euclidean case. A natural starting point could be Monma and Suris's paper [T7], where a 
partition of d-dimensional Euclidean space is constructed that has similar topology-limiting 
properties as the oriented Dirichlet cell partition. Another future research topic is an 
analysis of the practicality of implementing our algorithm for the Euclidean 1-bottleneck 
Steiner tree problem. We are also considering the merits of utilising this algorithm as an 
iterative heuristic for the Euclidean ^-bottleneck Steiner tree problem. Our interest in this 
topic stems from its application to the deployment of wireless sensor networks. 
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